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MSX Mid infra red imaging of massive star birth environments. II: 
Giant H ii regions 
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ABSTRACT 

We conduct a Galactic census of Giant Hll regions, based on the all sky 6cm dataset of 
Kuchar & Clark, plus the kinematic distances obtained by Russeil. From an inspection of 
mid-IR Mid-course Space Experiment (Af SX) and far-IR IRAS Sky Survey Atlas images 
we identify a total of 56 GH II regions in the Milky Way, of which 15% (65%) can be seen at 
optical (near-IR) wavelengths. The mid to far-IR fluxes from each GH II region are measured, 
and sample the thermal emission from the ubiquitous dust present within the exciting clusters 
of OB stars, arising from the integrated luminosity of the hot stars heating the cluster dust, for 
which we obtain log_L(IR)=5.5-7.3LQ. The mid-IR 21/xm spatial morphology is presented 
for each GH II region, and often indicates multiple emission sources, suggesting complicated 
cluster formation. IR colour-colour diagrams are presented, providing information about the 
temperature distribution and optical depth of the dust. For the clusters of our study, the dust 
is not optically thick to all the stellar radiation, thus the measured infra-red luminosity is 
lower than the Lboi- As the dust environment of a cluster begins to dissipate, the thermal 
emission and its optical depth ought to decrease even before the stars appreciably evolve. 
We see evidence of this in our empirical relationship between the integrated IR and Lyman 
continuum luminosities. 
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1 INTRODUCTION 

Massive stars, with O or early B spectral types on the main se- 
quence are, for the most part, born in giant molecular clouds 
(GMCs). Initial "clumps" made up of dense molecular gas along 
with appreciable concentrations of dust comprise the birth material. 
Once a luminous star becomes sufficiently hot, substantial Lyman 
continuum luminosity is produced and the surrounding hydrogen is 
rapidly ionized into a Stromgren Sphere. During these birth phases, 
the medium surrounding the stellar embryo is sufficiently dense 
10^cm~^)thattheHlI region radius is a small fraction of a par- 
sec, ~ 0.01-0.1 pc. This ultra-compact HlI (UCHll) region (e.g. 
Wood & Churchwell 1989) is detectable at cm wavelengths. It is 
surrounded by a natal dust "cocoon", which is somewhat opaque to 
(roughly) near infrared (NIR) and shorter wavelength radiation and 
so is heated by the star. These cocoons are observable as thermal 
sources in the mid infrared (MIR) and at longer wavelengths as has 
been shown previously (e.g. Crowther & Conti 2003 - hereafter Pa- 
per I). Further evolution of the individual stellar birth cocoons pro- 
ceeds as follows: The intense radiation of the star (likely aided by 
its wind) dissipates and evacuates the dust that gradually expands. 
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As it does so, its optical depth diminishes and the OB type excit- 
ing star becomes revealed. The UC II region also becomes larger, 
forming a compact H II and finally a normal H II region. 

For OB stars formed in a cluster, additional physics comes into 
play as the environment also contains gas and dust. Each luminous 
star affects this cluster material by heating the dust and gradually 
dissipating it by the strong stellar radiation fields and winds. In ad- 
dition, the Lyman continuum luminosity from a large ensemble of 
OB stars can ionize a substantial volume of space, resulting in an 
observable giant H II (GH II) region. In a few cases (e.g. W49A) in- 
dividual UCH II regions can contribute to the overall ionization of 
the cluster environment (e.g. Conti & Blum 2002, Alves & Home- 
ier 2003). We would expect that a GH II region would be ionized by 
OB stars which are, for the most part, already beyond the UCH II 
stage (e.g. Blum et al. 2002). We are interested in measuring param- 
eters of the dust content of the OB star clusters since that lifetime 
is much longer than that of the individual natal dust cocoons sur- 
rounding the OB star members, but may still be shorter than the 
stellar evolution time scale. 

We began this study by considering those Galactic (GH II) ob- 
jects that emit a Lyman continuum luminosity - N(LyC) - of at 
least IC" photon s"'^, thus excited by more than ten 07V* stars 
(O star equivalents - Vacca 1994). Such regions ought to repre- 
sent the main locations of massive stars in our Galaxy. Following 
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the technique of Smith et al. (1978) we re-evaluate the census of 
Galactic GH II regions based on modern radio surveys (e.g. Con- 
don et al.l989) and kinematical distances (Russell 2003), that is, 
they are derived from a Galactic rotation model and Hi line ve- 
locity measurements.^ This GHlI sample has formed the basis of 
an ongoing systematic NIR imaging and spectroscopic study of the 
stellar content of these objects (e.g. Blum 2003 and § 4). It is al- 
ready clear from our published and unpublished NIR imaging that 
many of the radio selected GH II sample contain identified OB star 
clusters. 

A primary objective of this paper is to study the newly selected 
GH II regions in the MIR with MSX at a higher (circular) spatial 
resolution than was possible with IRAS ^ . Kraemer et al. (2003) 
present studies of nearby star forming regions (Orion, Pleiades etc.) 
with MSX, whilst Spitzer has already began producing yet higher 
spatial resolution MIR maps (Churchwell et al. 2004; Whitney et 
al. 2004). Nevertheless, the present work extends previous studies 
to a complete census of the most luminous radio selected, mas- 
sive star forming regions. In many cases the spatial morphology is 
quite complicated, indicating the presence of more than one dust 
emission source, so previous H II region studies making use of the 
IRAS Point Source Catalogue (e.g. Codella et al. 1994) may be 
inappropriate due to the widely variable spatial scale of Galactic 
GHlI regions (Chan & Fich 1995). We will make measurements 
of the various IR fluxes through the IRAS filters using the Sky 
Survey Atlas (ISSA). 

We would like to compare measurements of the thermal emis- 
sion from the integrated cluster dust with the radio free-free emis- 
sion from the hot stars within. The former depends upon the far UV 
(FUV) luminosity of the ensemble of OB stars in the cluster(s) and 
the extent and temperature distribution of the dust within it. The 
latter depends upon the integrated extreme UV (EUV) luminosity 
of these same stars. In our Paper I we have proposed that these 
same measurements fit reasonably well a simple spherical model 
for the (mostly single star excited) UCH II regions. We will revisit 
this relationship below. 

It is well known that UCH II regions have spectral energy dis- 
tributions (SEDs) which are not fit by a single temperature black 
body (e.g. Wolfire & Churchwell 1994). Rather the observed SEDs 
of UCH II regions are broad with a peak near 100/im, correspond- 
ing to a dust temperature of 30 deg K. However, other hotter ma- 
terial is also present. What is the situation for the dust in those OB 
star clusters with dust? Here MIR and far-IR (FIR) photometry will 
help constrain the SEDs of dusty OB clusters. 

In Section |2| we present our revised GHlI catalogue, which 
is complemented by the largest GHlI regions in the Magellanic 
Clouds. The procedure for analysing the observations is contained 
in Section|3| We discuss in Section|4|the individual MIR images of 
this set of data which are presented in the Appendix. In Section|5| 
we discuss the similarities and differences among the objects us- 
ing different IR flux ratios (colours). We compare the GHll/HlI 
colours with those of the UCH II regions from Paper I. We also 
consider the relationships between the IR luminosity and the radio 
fluxes for all the data. A discussion is given in Section|6| 



The rotation model assumes axi-symmetric circular motions (no peculiai' 
velocities). Inside the solar circle a near/far side ambiguity is present which 
must be resolved with other data. 

^ The IRAS central wavelengths of 12/^m, 25 /^m, 60/xm and 100/^m have 
rectangular spatial resolutions of 30", 30", 1' and 2', respectively. 
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Figure 1. Distribution of GH II regions (filled circles) in the Galactic disk 
used in this study. The location of the Sun is indicated at a Galactocentric 
distance of 8kpc (Reid 1993), as ai'e radii of 4, 8 and 12 kpc (dotted lines). 



2 SELECTION OF SAMPLE 

In this paper we look at a planopy of massive star birth environ- 
ments in multiple excited GH II regions. At an initial stage of this 
study we based our sample on forty GH II regions of Blum (2003). 
This sample was directly drawn from Smith et al. (1978), which 
was based on low spatial resolution 6cm radio fluxes taken from 
Reifenstein et al. (1970) and Wilson et al. (1970), except that a 
revised Solar galactocentric distance of 8kpc was adopted (Reid 
1993). 

In preference to these older compilations, we used the higher 
spatial resolution 6cm study of H II regions by Kuchar & Clark 
(1997), based on the NRAO 9Im Condon et al. (1989) survey plus 
the ATNF 64m (Parkes) Condon et al. (1993) survey. These were 
used in preference to (admittedly yet higher) spatial resolution ra- 
dio data such as the NRAO VLA Sky Survey (NVSS) at 1.4GHz, 
since we required all sky coverage at an optically thin radio fre- 
quency, for which 5GHz is superior to 1.4GHz. In two cases (Ml 7 
and W43) we used 6cm radio fluxes from Downes et al. (1980), 
since these were missing from the Kuchar & Clark (1997) compi- 
lation, presumably due to saturation effects. 

From our initial list of candidate GHlI regions based solely 
upon Russeil (2003) radio derived kinematical distances, we first 
revised distances for several well-known regions based on more 
recent radio studies by Araya et al. (2002), Watson et al. (2003), 
Sewilo et al. (2004) and Corbel & Eikenberry (2004). This revi- 
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sion led several candidates to be revised to the near (rather than 
far) distance, demoting them to normal Hll regions, including 
G34.255+0.144 (Araya et al. 2002), GlO.617-0.384 and G10.664- 
0.467 (Corbel & Eikenberry 2004). 

Subsequently, further candidates were rejected on the basis 
of spatial location in the Milky Way inferred from the kinemati- 
cal distance assumed. Two were rejected since their Galactocentric 
distance was predicted to exceed 12kpc (based on a Solar Galac- 
tocentric distance of 8kpc, Reid 1993), namely G7. 472+0.060 and 
G173. 599+2. 798. In addition, two further candidates were omit- 
ted since their 'far' kinematical distance implied a distance of 
z~300pc from the Galactic plane, namely G351. 613-1. 270 and 
G351. 695-1. 165. It is likely that these two sources lie on the 'near' 
side of the Galactic Centre, such that they are normal H 11 regions. 
Two other candidates lying z ~200 pc from the Galactic plane are 
included (G70.300-1.600, G79.293+ 1.296) since they lie outside 
the Solar circle, and their distances are probably genuine. 

Finally, ten further candidates, some associated with the W49 
and W51 complexes, G2.901-0.006, G13.381+0.071, G26.103- 
0.069, G289.760-I.I55, G338.086+0.007, G342.382-0.044 were 
also dropped since their llfim fluxes were less than 5% of that for 
the mean GH II region, after being adjusted to a uniform distance 
in all cases. Finally, we also omitted the Galactic Centre cluster 
(G359.95 1-0.037) since we considered a study of its IR/radio prop- 
erties might be affected by unusual conditions relative to normal 
star forming regions. 

Table 1 gives the listing of the remaining 56 GH II regions as 
derived from radio measurements. Figure Q indicates the Galactic 
distribution of GH II regions in Galaxy, as viewed from above, pro- 
jected onto the disk, and along the disk of the Milky Way. This 
is analogous to Fig. 11 from Georgelin & Georgelin (1976), al- 
beit based solely on radio selected GH II regions. From the original 
Smith et al. (1978) sample, there was a lack of radio GH II regions 
on the far side of the Galaxy. It is apparent that more contemporary 
distance estimates has removed this puzzle. The spiral arm struc- 
ture is not immediately apparent (see Russeil 2003), although any 
interpretations must be treated with caution as the distances have 
been determined kinematically. The handful of objects so far stud- 
ied photometrically and spectroscopically in the NIR (e.g. see Blum 
et al. 2002) have spectroscopic parallaxes which move them closer 
to the Sun, two substantially. The appearance of this figure might 
change dramatically once more of the OB star clusters exciting the 
GH II have been identified and their distances determined by this 
method. With the exception of several sources at the Galactic Cen- 
tre itself, there is only one region apparently within a galactocentric 
distance of 4kpc, namely G347.61 1+0.204. 



preserved. A total ofl680 95.l'x95.l' images were created in each 
of the spectral bands, spaced 90' apart in each coordinate, provid- 
ing 10' overlap on adjacent images. Each image provides radiances 
on a grid in Galactic latitude and longitude, with a grid spacing of 
6". The image units are of in-band radiance (W m^"^ sr^^). 

Of the available bands. Band A is the most sensitive, but its 
filter footprint of 6.8-10. 8/im includes the astrophysically strong, 
broad silicate band at 9.7/^m and the 7.7, 8.6/im PAH features (e.g. 
Faison et al. 1998). We therefore selected the longer wavelength 
Band C and Band E images, which have isophotal central wave- 
lengths of Ac=12.13/im and 21.34/.im, respectively, with bandwidth 
of 1.72/im and 6.24/im. Consequently, MSX Bands C and E are 
narrower, higher spatial resolution analogues of the well known 
IRAS 12/im and 25/im bands. This is illustrated in Fig. 1 of Pa- 
per I. MSX Band D was also neglected, since its addition would 
not have added to the scientific results of this work. The calibration 
and photometric accuracy of MSX are discussed in detail by Egan 
et al. (1999). The zero magnitude flux is based on the Kurucz model 
for Vega (Cohen et al. 1992), and corresponds to 9.259x 10"" W 
m"^ (Band C) and 3.555 x lO"^'^ W m"^ (Band E). 

Individual GHlI and Hll region images were 
obtained from the MSX Image Server at IPAC 
http : 77irsa ■ ipac ■ caltech . edu/ applications/ 
MSX/. For our purposes, we sought integrated MIR fluxes 
in Janskys for comparison with radio fluxes. Consequently, the 
MSX point source catalogue was not appropriate, given the need 
to use variable apertures for these extended sources. Extracted 
MIR fluxes of each GHlI region are obtained from integrated 
fluxes using a circular aperture of diameter 6, as derived from 
MSX 2I/im images (see Table 2), with a mean diameter of 
6(21/im)=4 arcmin (40 pixels). Identical aperture diameters were 
adopted for the 12/im flux measurements. In all cases the apertures 
were centred on the radio positions, with the exception of a few 
spatially large nearby regions such as Ml 7, for which an aperture 
of 9= 9', centred at G15.062-0.693 was required to accommodate 
the entire GH II region. Thin annular rings were taken in all cases, 
in order to set the background level, taking care to avoid MIR 
sources in these annuli. 

Spatial integration was measured in GAIA (Draper, Gray & 
Berry 2001) using appropriate apertures, and provides fluxes in 
units of W m~^, after correction for the 6" x 6" pixel-area, i.e. 
8.4615x10"^" sr. Division by the filter bandwidth then provides 
fluxes in W m~^ Hz~^ (= 10^^ Jy). Finally, a multiplication by 
1.113 is required to convert the square arc pixels into the correct 
Gaussian area (Cohen 2002). The total scale factor for Band C cor- 
responds to a multiplicative factor of 27.45 in order to convert the 
IPAC integrated fluxes to mjy. 



3 MID AND FAR-IR OBSERVATIONS 
3.1 MSX Observations 

The Spatial Infrared Imaging Telescope (SPIRIT III) aboard the 
U.S. Department of Defense MSX satellite surveyed the entire 
Galactic plane in four MIR spectral bands, named A, C, D, and E, 
between 6 and 25/im at a (^3 pixel) spatial resolution of -^18. 3". 
A description of SPIRIT III is given by Price et al. (2001). The en- 
tire area within ±4.5° of the Galactic plane was surveyed at least 
twice, with four-fold coverage obtained in the first and fourth quad- 
rants, and to ±3° for the remaining quadrants. The redundancy was 
sufficient to permit combining the datasets onto a uniformly spaced 
grid, such that the inherent spatial resolution of the instrument was 



3.2 IRAS Observations 

Since the Hll regions under discussion here are extended, 
use of the IRAS Point Source Catalog (PSC) is generally 
not appropriate to our sample (e.g. as used by Codella et al. 
1994). Consequently, we extracted images for each Hll re- 
gion for the IRAS 25/^m, 60/im and lOO/im filters using 
the IRAS Sky Survey Atlas (ISSA) server, also at IPAC 
http://irsa. ipac. caltech. edu /applications/ 
ISSA/. In contrast with MSX, IRAS was not designed to 
achieve absolute surface brightnesses. Nevertheless, ISSA was op- 
timized to permit flux extraction for ('--^arc minute sized) extended 
sources (Gautier 2003). For comparison, Chan & Fich (1995) used 
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Table 1. Catalogue of Galactic GH II regions based on 6cm fluxes of Kuchar & Clark (1997) and distances from Russeil (2003), except where noted. Lyman 
continuum fluxes are corrected for dust extinction, following Smith et al. (1978) from which we also derive . We indicate whether the GH II region is seen 
optically (•) or in the NIR (o). 
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(1) Wilson et al. (1979); (2) Downes et al. (1980); (3) Israel (1980); (4) CasweU & Haynes (1987); (5) Kuchar & Clark (1997); (6) Araya et al. (2002); (7) 
Russeil (2003); (8) Watson et al. (2003); (9) Sewilo et al. (2004); (10) Corbel & Eikenberry (2004); (11) fliough see Blum & Damineli (1999); (12) Madore & 

Freedman (1998); (13) Harries et al. (2003); (14) Wink et al. (1983) 
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earlier Full REsolution Survey CO-adder (FRESCO) data products 
for their study of H ii regions in the range 70< I < 235. 

Again, GAIA was used to spatially integrate the ISSA images. 
The apertures were generally centred on the radio peak with aper- 
tures selected to measure the entire flux associated with the GH II 
region in the 60/im filter, whilst attempting to neglect the contribu- 
tion from neighbouring regions. Identical apertures were adopted 
for other filters. Since ISSA pixels are uniformly 1 .5 arcmin across, 
ISSA apertures were correspondingly larger than MSX aper- 
tures, with a mean diameter of ^(60/im)=10.5 arcmin (7 pixels). 
ISSA images provide fluxes in MJysr"'^. The solid angle of one 
ISSA pixel is (1.5/(57.3 x 60))^ = 1.90386 xlO"'^, such that one 
needs merely to multiply spatially summed counts by a factor of 
0. 190386 to obtain integrated fluxes in Jy. We have compared ISSA 
25^im to MSX 21/im fluxes in all cases, revealing an average ratio 
of 1.7±0.8. One would expect the IRAS fluxes to be greater owing 
to the greater bandpass and lower spatial resolution. 



4 MIRMOPHOLOGIES 

In the Appendix, available electronically. Figs. Al-10 of MSX 
Band E (21/im) images are presented of all Galactic and Magel- 
lanic Cloud GH II regions. Each of these figures contains six ob- 
jects, with a field of view of 15 x 15 arcmin. The intensity is in log- 
arithmic units (W m~^ sr"'^), scaled between the brightest knots 
and those 1.0 dex fainter, as indicated in the adjacent key. In these 
figures one is looking at the integrated dust emission, which comes 
from unresolved sources, resolved objects (the cluster dust), and 
from an extended halo associated with the formation process. In Ta- 
ble 2, we indicate whether individual sources are single or multiple, 
compact or elongated, or are associated with diffuse components. 

All of the sources are resolved with MSX, that is, their 
sizes are greater than 18", the spatial resolution of the instrument. 
This corresponds to different physical scales, depending on the 
source distances, as may be seen in the images. Many of them 
have multiple components, suggesting several star forming clusters 
are present within the ionized hydrogen regions (e.g. W49A). The 
sources in this paper frequently show diffuse emission surrounding 
them, often with loops and tails. This is in contrast to the single 
star excited UCH ii regions of Paper I where such features were 
less common. Unattributed comments concerning NIR photome- 
try of individual sources are from unpublished observations (Blum, 
private communication). In many cases, the MIR MSX images 
closely mimic radio continuum datasets, also common to UCH ii 
regions from Paper I, and the study of Cohen & Green (2001). 

We now discuss a subset of the GH ii regions, grouped accord- 
ing to whether they are observed optically or in the NIR. 

4.1 Optically visible GH 1 1 regions 

Within the Milky Way, ~15% of the GHii regions are optically 
visible - namely NGC 3372 (Carina), M8, M17 (Omega), RCW42, 
RCW49, RCW97, NGC 3576 and NGC 3603. Their MIR mor- 
phologies are complex, ranging from a peculiar V-shape of the dust 
in the case of M17, which contains a luminous cluster (Hanson et 
al. 1997), to several distinct regions of very bright dust emission in 
Carina (Rathbome et al. 2004) and NGC 3603 (Numberger 2002). 

Rathborne et al (2004) present a detailed morphological dis- 
cussion of Carina, such that it is not discussed here, except that the 
measurements presented here are centred on the Northern Molec- 
ular Cloud, containing the cluster Tr 14, such that rj Carina is 



excluded. For NGC 3603, both the brighter, compact north-west 
source and extended SW source are included in our MSX mea- 
surements. The optically visible cluster lies between these two IR 
sources slightly to the north, the stellar content of which is well 
studied (e.g. Moffat 1983), and contains three H-rich Wolf-Rayet 
stars (Drissen et al. 1995). Finally, a NIR image of the star cluster 
associated with NGC 3576 is given by Figueredo et al. (2002), who 
also include K band spectroscopy. 

A strong elongated MIR source is found at the radio position 
of RCW49, with a similar morphology to M17. A cluster is present 
in NIR images, whilst the Westerlund 2 cluster lies immediately 
to the north of the peak MIR emission, which hosts two WR stars 
(Shara et al. 1991). Finally, RCW42 reveals a bright MIR source, 
with an extended halo, which is spatially coincident with a NIR 
cluster. 

Turning to extragalactic sources, 30 Doradus, alias NGC 2070 
(Tarantula Nebula), is the prototype optically visible GHii region 
of the Local Group. Spectroscopy of stars within the central clus- 
ter, R136, reveals a large concentration of early O stars (Massey & 
Hunter 1998), which also hosts a number of Wolf-Rayet stars (e.g. 
Walbom 1973; Crowther & Dessart 1998). In the MSX images, 
dust emission appears as an incomplete torus with an opening to the 
east. NGC 346 is rather less massive, although it hosts many young 
OB stars (Massey et al. 1989). Although fairly weak, NGC 346 is 
spatially extended NW to SE in the MSX image. 



4.2 NIR GHii regions 

Approximately 2/3 of the GHlI regions can be readily identified 
on K-band 2MASS images via diffuse nebulosity. AISX images of 
W3 have already been presented by Kraemer et al. (2003), and so is 
not discussed further here. In addition, several sources host UCH II 
regions, including G29.944-0.042, G70.300-i-l.600 (K3-50A) and 
G79.293-f 1.296 (DR7), as discussed in Paper I, the latter hosting a 
NIR cluster (Dutra & Bica 2001). 

GlO.159-^.349 (W31), W42 and G333.610-0.217 have a sim- 
ilar MIR spatial appearance, with a compact spherical core sur- 
rounded by an extended halo. Each host a luminous cluster as re- 
vealed by NIR images of Blum et al. (2000, 2001, 2003). In con- 
trast, W43 includes three separate dust emission sources, of which 
the central source contains a young cluster, itself hosting a Wolf- 
Rayet star at its core. W49A is more complicate still, revealing 
several distinct distinct dust emission sources that correspond to in- 
dividual clusters of OB stars (Conti & Blum 2002). Note the MIR 
spatial morphology of W49A is exactly matched by Br7 emission 
on the deep NIR image of this field by Alves & Homeier (2003). 
GH II regions within the W51 complex exhibit similar spatial mor- 
phologies. 

G298.227-0.340 and G35 1.467-0.462 have a rather simple 
MIR morphology, with compact cores lacking an extended halo. 
These source, alias IRAS 12073-6233 and IRAS 17221-3619, 
were included in the IR/radio study of Hll regions by Martin- 
Hernandez et al. (2003). Finally, the compact NIR H ii region iden- 
tified in Gl. 149-0.062 by Blum & DamineU (1999) is spatially co- 
incident with the dominant MIR source Sgr D at Gl. 128-0. 104, 3' 
SE from the radio peak. 



4.3 MIR GHii regions 

The remanding ~l/3 of the GHii region sample exhibit neither 
optical nor NIR signatures associated with the radio/MIR source. 
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Table 2. MSX 12//m (Band C) and (Band E) and IRAS ISSA (60/im and lOO/xm) fluxes for radio selected Galactic GH II regions, including a brief 
summary of the 21/im spatial morphology, namely single (s) vs multiple (m), compact (c) vs elongated (e), containing diffuse emission (d) or not. 
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For this subsample, the extinction to the ionizing source is ei- 
ther rather high or in a few cases identification as GH ii regions 
may be in question. MIR morphologies for this group span a wide 
range. In some cases the GH ii regions reveal a compact dust source 
with associated extended halo (e.g. G345. 555-0.043) that resem- 
ble the morphology of regions with luminous NIR clusters such as 
W31. In other cases bright, extended dust sources are observed (e.g. 
G2.303-i-0.243), reminiscent of RCW49. Several exhibit morpholo- 
gies involving multiple point sources plus diffuse emission (e.g. 
W24, G20.733-0.087). Finally, several show primarily diffuse dust 
emission (e.g. GO.489-0.668, G3.270-0.101). 



4.4 Colour-colour diagrams 

The dust emission of the GH II regions arises primarily from the in- 
terstellar material within or close to the exciting star clusters. Some 
emission may also have an origin from individual circumstellar dust 
cocoons surrounding UCH ii regions (or hot cores) within the clus- 
ters. As with stellar magnitudes UBV or JHK, it will be instructive 
to plot analogous colors for the MIR or FIR fluxes. In Figure 2(a) 
we make IR colour-colour plots using the log flux ratio (colour) 
100/umy60/^m versus 60/im/25/xm for the UCH ii of Paper I (open) 
and for the GH ii regions of this paper (filled). Similar plots for a 
smaller number of normal H II regions were presented by Chan & 
Fich(1995). 

We see immediately that the sources in the two plots have very 
similar colours, although there is more scatter in the UCH ii points. 
Recall that Wood & Churchwell (1989) used colour-colour plots 
similar to these to identify UCH II regions. We see that such criteria 
will also pick out GH ll regions. The outliers in Fig 2(a) with the 
lowest lOO/im/60/xm ratios are GO.361-0.780 and GO.394-0.540. 

To get a feeling for the dust temperatures associated with these 
colours, we used the black body distributions convolved with the 
IRAS filter profiles (based on Beichman et al. 1988). In each plot 
there is a main body of points with outliers, the most extreme 
of which are usually those with complicated morphology. In both 
cases we see that on average the 60/25/im colours lead to higher 
temperatures than the 100/60/im colours. This is a quantitative way 
to demonstrate that the SED of the dust of UCH II, GH II and H Il's 
do not follow simple black body laws. Note that the main body of 
points of the UCH ii are at slightly cooler temperatures but there is 
plenty of overlap. 

Using the MSX data we have also made an analogous plot of 
21/12/Ltm versus the IRAS 60/25/im colors, which is presented in 
Fig. 2(b). The inferred temperatures were still higher with these 
shorter wavelength MSX colours, although we acknowledge that 
the 12/im fluxes are influenced by the silicate absorption feature 
which affects this wavelength. The observed scatter for the complex 
morphology UCH ii regions (open triangles) was much greater than 
for the simple UCHii regions (open squares) and GHii regions 
(filled circles). 



5 RELATION BETWEEN RADIO EMISSION AND FIR 
FLUXES 

The radio emission of H ii regions arises from the integrated Ly- 
man continuum luminosity (EUV) of the hot star content of the 
cluster(s) within. What is measured is the free-free emission of the 
ionized hydrogen and/or the strength of high level recombination 
hydrogen emission lines, typically obtained at cm wavelengths. 
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Figure 2. Colour-colour diagrams: (a): IRAS log flux ratio (colour) 
1 ()0/im/60^tm versus 60fim/25fim for the GHll regions newly presented 
here (filled circles) and UCH II regions from Paper I (open squares: sim- 
ple morphology, open triangles: complex morphology), (b): log flux ratio 
(colour) of IRAS 60/um/25;im versus MSX 21/im/12/im for the GHll 
regions (filled circles). 



These can be converted to the number of Lyman continuum pho- 
tons (NLyC) producing the emission using well known physical 
relationships, if the distance is known. What one wishes to obtain 
is the NLyC photons emitted by the stars, so that one may infer 
their stellar properties. The observed NLyC could be affected by 
leakage of photons from the region without ionizing the material, a 
so-called density bounded nebulosity. Some emitted photons might 
also be absorbed by the dust within the nebulae. In both cases, the 
observed NLyC would be smaller than the emitted NLyC and one 
would underestimate the hot star content. How important are these 
effects for GH ii regions? 

UCH II regions have very dense cocoons overlaying the ion- 
ized hydrogen but also tend to have extended hydrogen emission 
halos, suggesting that their inner environments are density bounded 
and/or clumping may play a role (Ignace & Churchwell 2003) and 
LyC photons are escaping. Whether clusters are density bounded is 
a longstanding problem that is not fully resolved, but does not need 
to be settled here. It would be important to resolve this issue were 
we to try to infer the stellar content of our clusters from the coimted 
NLyC photons alone. 
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We now turn to the role of the dust. Its emission is related to 
the quantity and quality of material present, and the heating by the 
hot stars in the clusters. Their emitted radiation peaks in the far 
UV beyond the Lyman limit; the thermal dust emission peaks near 
100/im. We have estimated the total IR integrated luminosity (Lir) 
of the dust using a standard relation (e.g. Chan & Fich 1995), with a 
minor modification to include the more reUable MSX MIR fluxes. 
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100 )\ 



If the dust is optically thick to stellar radiation, the Lir is a good 
estimate of the hboi- With optically thin dust, stellar radiation es- 
capes to space and is not available to heat the material and the LIR 
will underestimate the total luminosity Lboi. 

In Fig. 3, we plot the observed Lir inferred from MSX and 
ISSA, versus log N(LyC), which are both distant dependent. This 
ratio is distance independent but strongly spectral type dependent 
(e.g., Churchwell 2002). (In Paper I, we made an analogous plot 
but used the lOO/nm limiinosity for the ordinate. That relationship 
is similar to that shown here.) The UCHii are to the right at lower 
log N(LyC); the GH II to the left. H II regions would fill the space 
immediately to the right of the GHii. The lines represent model 
predictions from Lumsden et al. (2003) for clusters with a nominal 
Salpeter IMF (solid line: upper) and single stars (dotted line: lower) 
based on CoStar stellar models of Schaerer et al. (1996). Similar 
results regarding Lyman continuum ionizing photons are obtained 
with alternative WM-basic models (Smith et al. 2002). Whilst the 
majority of the UCH ii regions lie between the cluster and single 
star predictions, the GH II regions actually lie below the (extrapo- 
lated) single star predictions - their IR luminosity is a factor of a 
few too low versus the cluster predictions. What is going on here? 
We assert that the problem is that the models assume all the stellar 
radiation is absorbed by the dust, which is optically thick, thus Lir 
= Lfcoi. It could be that some of the UCH II regions, which lie close 
to, or just above, their model predictions, are optically thick. But 
the clusters are not. 

What are we to make of the large nimiber of UCH ii regions 
between the single star and cluster predictions? The most reason- 
able explanation is that these objects contain extra stars which con- 
tribute to the Lir but not so much to the N(LyC) numbers. They 
behave, therefore, as miniclusters. This explanation was also given 
by Lumsden et al. (2003). But there are still a dozen or so UCH ii 
above the cluster predictions. This cannot be an optical depth ef- 
fect as smaller values would push the Lir lower. It is possible that 
some of these sources have either dust absorption of the LyC pho- 
tons and/or they are density bounded thus they would be moved to 
the left. Indeed, even some of the UCH II regions located between 
the model predictions could have these problems. Note that a cou- 
ple of UCH II regions are below the single star predictions. These 
might not be optically thick. 

Let us now turn to understanding the placement of the GH ii 
regions in Fig 3. There is, curiously, less scatter in the vertical dis- 
persion than for the UCH II regions (one would have expected more 
scatter given the larger number of free parameters for the clusters). 
We can assert that those clusters along the top edge of the rela- 
tionship are on average dustier, or have a higher optical depth, than 
those along the lower boundary. The vertical or horizontal scatter is 
a factor of a few. While we may be miscounting the N(LyC) emit- 
ted by the stars which would move them to the left, we suspect 
that the optical depth effect in the dust is the more important vari- 
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Figure 3. Comparison between Lyman continuum log flux inferred from 
radio free-free emission and IR luminosity of UCH II (from Paper I, open 
symbols) plus Galactic and Magellanic Cloud GHUii regions (filled sym- 
bols) inferred from ISSA and MSX. The solid Unes are the predictions 
of the CoStar models from Lumsden et al. (2003) for clusters, whilst dotted 
line are for single stars. 



able here.^ Included among the lower boundary group is NGC 346, 
which contains evolved O type stars, thus is somewhat older than 
the others. The vertical range of the clusters of Fig. 3 at a given 
N(LyC) could then be mostly an age effect as the dust is gradually 
being dissipated and the optical depth decreasing. The upper edge 
of the cluster distribution is a factor of 2 or so below the model 
predictions. This gives an estimate of the overall dust optical depth 
to stellar radiation, which would be of order r=0.7. 



6 SUMMARY 

We have carried out a complete census of 6cm selected H ii regions, 
for which 56 giant H II regions are identified, based in part on mid 
and far-IR fluxes from MSX and IRAS. Of course, the radio se- 
lected GH II catalogue presented here is not complete. In addition, 
there are several clusters close to the centre of the Galaxy with 
ionizing fluxes in excess of 10"'" s^^. These include the Arches 
(Serabyn et al. 1998) and Quintuplet (Figer et al. 1999) clusters. In 
addition, older massive clusters that have already lost the bulk of 
their associated gas are also omitted (e.g. Westerlvmd 1: Clark & 
Negueruela 2002). 

We present 21/im images of all giant Hll regions using 
MSX. This wavelength is sampling the thermal emission from the 
dust in the vicinity of these OB clusters, material left over from the 
birth of these objects. The spatial morphology of most of the dust 
emission is complex; many sources are double or multiple, suggest- 
ing several common sites of recent star formation. Extended diffuse 
dust emission is also common, in contrast to our analogous study 
of UCH II regions (Paper I) where such features were rare. In the 
few cases in which an UCH ii region is found within a larger GH ii 
region, the dust from the former entity dominates the emission; it 
is only a small part of the measured NLyC. 

We have also measured MIR and FIR emission from the 



Errors in the distance affect both axes similarly. 
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ISSA archive of these same objects. We used the MSX sateUite 
data to interpret the FIR IRAS fluxes which were obtained at a 
lower spatial resolution, particularly the longest wavelength one at 
100/im. These wavelengths sample an appreciable fraction of the 
dust emission and can be used to infer properties of its SED and 
total luminosity. Tables 1 and 2 provide detailed information about 
each source, with relevant flux measurements. A careful examina- 
tion of the Tables indicates that 75% of the sources show a 100/xm 
flux larger than that at 60/im, thus the peak of the dust emission 
SED for the cluster excited regions is near the former wavelength 
or beyond. This was already known to be the case for UCH II re- 
gions but is a new result for the GH II regions. 

In the colour-colour diagrams we plot the MSX and IRAS 
flux ratios for the UCH II regions of Paper I and the luminous 
GHlI of this paper. There is a main body of points in each plot, 
with outliers, and while the UCH II are a bit redder there is con- 
siderable overlap in the colours with the cluster objects. Under an 
over simplified assumption that each colour can be represented by a 
black body temperature, we also show these numbers on the plots. 
The temperatures inferred from the colour are mostly lower for the 
longer wavelength colours in both cases. IR colour-colour selection 
criteria for UCH II regions will also pick up GH II and H II sources. 

We infer the FIR luminosity by weighting and combining the 
IR fluxes using the distances in the Tables; the N(LyC) is also listed 
there. In Fig. 3 we show the relationship between these two mea- 
sures. With some scatter we see that as the number of OB stars (the 
abscissa) increases, the dust emission (ordinate) does also (more 
stars). The UCH II points lie, for the most part, between the pre- 
dictions for single stars and those for small clusters. Some objects 
lie well above the cluster relationship; these most likely have some 
absorption or escape of LyC photons from the nebula (quite possi- 
ble if there is significant clumping (Ignace & Churchwell 2004). In 
both cases the points would move to the left towards the predicted 
line, sometimes substantially. 
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Figure Al. MSX Band E images of GH II regions. Each field in this and subsequent figures covers a field-of-view of 15 x 15 arcmin, and is presented in a 
logarithmic intensity scale (Wm~^ sr~^). The centre of each field refers to the quoted coordinates of the radio peak. 
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Figure A2. MSX Band E images of GH II regions. 
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Figure A3. MSX Band E images of GH II regions 
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Figure A4. MSX Band E images of GH II regions 
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Figure A5. MSX Band E images of GH II regions 
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Figure A6. MSX Band E images of GH II regions 
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Figure A7. MSX Band E images of GH II regions 
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Figure A8. MSX Band E images of GH II regions 
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Figure A9. MSX Band E images of GH II regions 
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Figure AlO. MSX Band E images of Galactic and Local Group GH ll regions. 
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